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An eflective-medium theory is proposed for random weakly nonlinear dielectric media. It is 
based on a new gaussian approximation for the probability distributions of the electric field in each 
component of a multi-phase composite. These distributions are computed to linear order from a 
Bruggeman-like self-consistent formula. The resulting effective-medium formula for the nonlinear 
medium reduces to Bruggeman's in the linear case. It is exact up to second order in a weak- 
disorder expansion, and close to the exact result in the dilute limit (in particular, it is exact for 
d = 1 and d = oo). In a high contrast situation, the noise exponents are k = k' = near the 
percolation threshold. Numerical results are provided for different weak nonlinearities. The use of 
the Bruggeman formula as a starting point for nonlinear homogeneization theories in dimensions 
d > 2 is questioned on the basis of known exact bounds on the noise exponents. 



PACS numbers: 77.84.Lf, 77.84.-s, 77.20. Ht, 05.40.-a 



I. INTRODUCTION 

Numerous papers have been-de voted to the problem of determining the effective constitutive law of random nonlinear 
composites in electrostaticajlij, of importance for understanding various nonlinear phenomena, either entering the 
class of "weakly nonlinear" (WNL) phenomena, £« that of "strongly nonlinear" phenomena such as random fuse-type 
or dielectric breakdown-type failure of materialsEZI. . .. . 

A convenient approach to the homogeneization problem is the "energetic" oneE^EJ, which can be summarized as 
follows. Let Wx{E) be the local potential (the energy density) from which derives the local nonlinear constitutive law. 
The homogeneization step consists in computing the effective potential W as the volume average 

W={w.,{E)), (1) 

from which the effective constitutive law can be deduced (see below). Denoting by Pa the volume proportion of 
component ("phase") a in a composite medium where Wx — Wa if x is in phase a, and introducing the volume 
averages per phase {■)a, W can be written 

W^ = ^p„(u;„(S))„. (2) 

a 

A class of effective-medium theories, mostly concerned with weak power-law nonliaearities, use our ability to 
compute the secpad moment of the electric field E in the various phases of the compositell3, completed by a decoupling 



approximatioratJ 



{w^iE))^^Wo.{{Eyj^). (3) 



This approximation is exact only when Wa is quadratic (the linear theory) , or when the field is uniform, which occurs 
only if all the Wa are equal, i.e. in a homogeneous medium. 

The purpose of this paper is to present a means to overcome this approximation, starting from the following remark: 
actually, with an additional hypothesis of ergodicity per phase, a volume phase average can be reinterpreted as an 
average over some probability distribution for the electric field in the phase under consideration. In this paper, volume 
averages will henceforth be identified to statistical averages. Then, the problem of computing W, Equ. reduces 
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to computing the probability distribution Va of the field in each phase, by means of which the phase averages (■)q 
can be carried out, without appealing to (^. 

Note that the full probability distribution 'P(E) of the electric field in the medium is then 

P(E)=^p„7'„(E). (4) 

Thus, the probability distributions of the field per phase are functions of great importance, since they cajxy all the 
information needed to build up the effective nonlinear ptpperties. A few analyses of 'P(E) are availablecj, mainly 
focused on its tail of relevance for breakdown phenomenaE3. However, up to my knowledge, none has been specifically 
aimed at the individjual components Va (E) . Chen and Sheng studied numerically the distribution of the local field 
in binary compositescJ. An analytical model of the same distrittution was subsequently proposed by Barthelemy and 
Orland, with excellent agreement with the latter simulationscj. Both works were limited to linear composites, in 
a weak-contrast situation. jA, first numerical attempt in the continuum, in a strong-contrast situation this time, is 
due to Cheng and TorquatcE3. They carried out two-dimensional numerical calculations on disordered systems with 
inclusions of various shapes. An exact calculation Im, V va a linear medium with Hashin-Shtrikman (HS) geometry 
was recently made using a density of states approachEa. In these works, the distribution V was found to be essentially 
bimodal for a binary medium (save in the HS|-seometry where fine-structure peaks are not wiped out by positional 
disorder), especially at high dielectric contrasted. This is understandable in wiew of the decomposition (0). 

Still, in parallel with the second moment of the field, the first moment in each phase (the average field) can be 
computed as well from any linear effective- medium theory, without requiring any a priori knowledge about Ep. Widely 
used in mechanics of random continuous media in the framework of the so-called "thermoelastic problem" E3, this fact 
seems to havej^caped most of the litterature on random dielectric media, excepted in the work by Ponte Castaneda 
and KailasamEj where it is used to compute reference fields around which non-linear potentials are expanded prior to 
homogeneization. Knowing the first two moments is enough to attempt a modelization of Va- The proposal examined 
in the present paper is to approximate each Va (E) by a gaussian vector distribution, and thereby, to estimate (^ for 
any potential, without appealing to the strong (and somewhat artificial) decoupling assumption (^). 

The scope of the flgper is limited however to WNL media, for which it is sufficient to compute the field distribution 
to linear order onlya'EI. An extension of these considerations to strongly nonlinear media will be examined elsewhere. 



II. GENERAL FRAMEWORK 



Let us therefore consider a d-dimensional random medium described by WNL thermodynamic potentials of the 
type 

w.{E)=wf\E)+wf{E), (5) 

where 

u}^{E) = \eE^ (6) 

{£{x) being the local permittivity entering the linear part), and w^{E) is the nonlinear (i.e. non-quadratic) part of 
the potential, assumed to be small. The local constitutive law is 

The potential is prescribed in each phase: w^iE) — Wa{E) if x £ a. Under the ergodic hypothesis introduced 
above, and a "boundary" conditioned (E) = E°, the effective potential reads 



where 



W{E") = Wun{E") + Wnl{E°), (8) 

W'ii„(£;") = (9) 

a 

w-„i(ii;") = ^p„«'(i?))„. (10) 
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The effective constitutive law in the homgeneized medium is then 



(11) 



To first order in tu"', the field distributions entering (||) and ( pO| ) need only be computed at the linear leveltla. The 
procedure is the following one. An extended linear cadsotropic homogeneization problem is first considered, with an 
artificial ferroelectric-like part added to the potentially: 



(12) 



where (x) = e^- and Pi{x) — P" if a: G a. The corresponding local constitutive law is Di = SijEj + Pi. Assuming 
one knows an estimate for 

WUE") = (*f(i?)), 

(cf . below) , the following exact relations hold for the distribution Va of the original linear problem 

1 dWnn 



ij —^^ij 
Pi=0 



(13) 

(14a) 
(14b) 



I emphasize that the derivatives are computed at the point Pi — 0, efj — e°'5ij where e" is the scalar permittivity in 
(|) and (|). There, w^™ = w^™ and Wiin = Wiin- Note that (|l4b[) is a slight generalization of the formulatSI 



2 dWy^ 



(15) 



Both (14a) and ( |l4q ) are particular cases of a more general theorem, the demonstration of which can be found in 
appendix B of the review by Ponte Castaheda and Suquet on nonlinear compositeaHj. This theorem is a consequence 
of the equality {D-iEi) — {Di){Ei) which holds only for particular boundary conditions, or in the infinite-volume limit 
if there exists an effective permittivity. 
Let us set 



(16) 
(17) 



with = M"£;0 and = C^E^E° + Cl{S.,j-E°E°), since the averages only depend on the direction E^ Ef/E° 
of the macroscopic field. I introduce now the gaussia-ii approximation for the probabiHty distributions l^a 



[(27r)''det(C"; 



ll/2' 



(18) 



which constitutes the main ingredient of the theory. 

There remains to compute an approximation for the quantities {Ei)^ and {EiEj)^, in oprder to completely define 
the above gaussian distributions. Let us therefore estimate W\in{E'^) through an extension of the Bruggeman self- 
consistent effective- medium approximationEHl to potentials of the form (p^). The formula is given and discussed in 
Appendix A. Using the scalar permittivities and Pi — 0, this estimate reduces to 



WUE°) = -eoE°' 



(19) 



where Eq is the usual Bruggeman effective permittivityE2l. In addition, the derivatives (14a) and (14b) of the estimate 
for WiiniE^) yield 
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M" = /i"£;°, (20) 

i2 



(21) 



where 



- .^ + (ll)sO ' (2^) 
y=7^' (23) 



d + 2 



(24) 



Note that the Bruggeman equation for Eq reads (/i) = 1, or ahernatively {jie) — (/i)eo = So- From these relations, 
we check that: (i) {eE) — e^{E); (ii) (eE'^) = e^{E)'^, where (E) = E° , as must be. Expressions for and CJ_ are 
readily obtained from (|2l]). 

In Fig. 1 are displayed curves for the probability density V\\{E\\/ E^) of the normalized longitudinal electric field, 
= E • E", for a binary medium with dielectric ratio eijex = 5, for various fractions p2 of material 2. The space 
dimension is d = 2. With V{E) given by (§), and E written as E = £'||E" + Ej^ (which defines E^), V\\ is obtained as 



V\\{E\\IE^)^E^ Id'^-^E^Vm^E^Y.Vc. io (25) 

where Ex_ is the transverse component of the field: E^ = E — E||. The probability density tends towards a Dirac peak, 
centered at i?||/£'° = 1 as p2 — >■ or p2 — > 1, since the field is then equal to the applied field E°. For intermediate 
volume fractions, the bimodal character of the probability distribution is evident. 

Under the same conditions, Fig. 2 displays the probability density of the scaled modulus of the transverse electric 
field, E^jE'^ . This distribution is defined, for E\_ > 0, by 

VAEJE0) = E0'-'J rii^||^»E.7^(E)^i^"^p. ^^^^^^-i^,^, (26) 

where Sd = /T{d/2) is the surface of a d-dimensional unit sphere. Unlike the previous one, this distribution 

is not multimodal since the mean value of the transverse field is (Ej^)q, = in each phase, whatever the volume 
fractions. 

Fig. 3 displays (i?||/-E") for a higher dielectric ratio £2/^1 = 1000. The concentrations are P2 = 0.12, 0.24, 0.36. 
These plots can be compared to that obtained from hard disks simulations by Cheng and Torquato (CT) in Fig. 3 of 
Ref.Ej, for the same dielectric contrast and concentrations p2 — 0.2, 0.4, 0.6. These concentrations differ from ours 
by a factor 0.6. Apart from these differences, the overall features (shape and heights) of the probability distributions 
are well rendered: a Dirac peak close to E\\ — indicates that the field is almost null in phase 2, and a widely 
spread contribution from phase 1 for the highest p2 indicates an enhancement of the fluctuations in this phase as the 
percolation threshold (for the Bruggeman theory), or the jamming threshold (for the CT simulations) is approached. 
The differences in the values for p2 can be explained by the fact that the present theory relies on the Bruggeman 
effective-medium formula, of relevance for cell-materials but not really adequate for hard disks. Moreover, the overall 
agreement between our distributions and that of CT is not so good at moderate dielectric contrast: Bruggeman's 
theory completed by the gaussian approximation somewhat overestimates, because of its percolating nature, the 
width of the probability distributions. Better adequation with the CT simulations would probably be obtained with 
an effective-medium theory for dielectric-coated inclusions, which would prevent percolation. Since our main objective 
however is to discuss an effective-medium therory for WNL composites, such an improvement will not be considered 
here. 

With these results in hand, analytical investigations of the theory can be carried out in some limiting cases, as well 
as numerical ones in more complicated situations. The next section examines particular WNL potentials. 
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III. WEAKLY NONLINEAR POWER-LAW POTENTIAL 



As a first application, a WNL local potential of the form 

w,{E)^]^eE^ + -^xE^ + 0{E''), (27) 

is considered, where the E'^ term is the first in a series of corrections to the linear behaviour. The corresponding 
constitutive law is 

D, = [e + xE"" + 0{E^)]E,. (28) 
This approximation actually is a weak-field one. 



A. Effective-medium formula 



The effective potential takes the form 



where is the Bruggeman result, ani 



i' = {xE')/E''^Y.Po^X.{E')jE'' 



(29) 



(30) 



The phase averages {E^)a are carried out with the help of Wick's theorem which allows to write down expressions 
for the integer moments of a centered vector gaussian distribution by mere inspections: with Ai = Ei — M", 



{AiAiAjAj) a — {AiAi)a{AjAj)a 2{AiAj)a{AiAj)c 



(31) 



We deduce 



{E'*')a — {EiEiEjEj)a 



+ (d- 1)C2 + Af"^ 



a2 



A few simplifications yield the simple result 

X" - [(1 + 2/d)y' + (2 - 2/d)z^ - 2] 



(32) 



(33) 



As a first check, we remark that in a non-disordered situation where e and x ^re constant in the medium, we have 
£a — e and /i = 1, so that y ^ z ^ 1; whence xo = Xj ^ '^^^ expected. 



B. Weak-contrast expansion 

A further check for the self-consistent formula ( ^3| ) consists in examining its weak-contrast limit. An exa«tjeptcpession 
is known for any nonlinear potential, which is first briefly reminded. In the weak-contrast expansionB'Lrll3ll3, the 
local potentials Wx{y) are assumed to fluctuate weakly around their mean value w''''(y) = {wx{y))- Introducing a 
bookkeeping parameter t to be set to 1 in the final results, the contrast, w'i'\y), is defined by 

u;,(2/) = u;(°)(2/) + u;«(2/)i, (34) 

so that {wjx'\y)) = 0. An expansion for the effective potential is sought for in the form: 

WiE'^)^W^°^ +w'^^h + w^^h'' + ■■■ (35) 
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Then! 



g 



W{E') = {w{E')) - ( ^)t' + 0{t% 



(36) 



where 



0\2 



S'd l + (r-l)(k-E0)2^ 

r = e||/£±. 



(37) 

(38) 
(39) 

This resuh is the exact one0. An analogous expapaion for the complementary potential W{D^) defined from the dual 
homogeneization problem can be written downllj'tSI. 

Let us therefore set e = (e) + Set and x = (x) + '^X t- The Bruggeman permittivity is exact to second order in 
the constrast: 



e° = is) 

Expanding ( |33| ) to second order in t entailfEi 



2(d + 8)(fc2) A{SeSx) 



did + 2) d (s){x) 



(40) 



(41) 



which can directly be obtained from a first-order expansion of ( p6| ) in x- 
Comparing now formula (^sh to the widely used approximation 



x'^X°2 = J2^a{ErjE''^y'{f.\), 



we see that the latter is not exact to second order in the contrast, save in = 1: its expansion indeed reads 

'2{d + 2){5e^) 4(&5x>' 



d2 d {e){x) 



t^'+Oit^) 



(42) 



(43) 



Hence, the gaussian decoupling which accounts for the vector character of the electric field, of importance in nonlinear 
problems, is superior to the simple approximation {E'^) ~ (i?^)^. 



C. Dilute limit 

For a binary medium with two components of constitutive parameters (£i,Xi) ^-nd (e2,X2), the dilute limit is the 
limiting situation where the volume fraction p of (e.g.) component 2 is small. Setting 



£2 - £1 



an expansion, to first order in p of (|33|) yields 

x"-Xi + ((X2-Xi)(T-l)^ + XiT2 



£2 + (rf - l)£l 

d(rf + 8) 



d + 2 



p + 0(p^). 



(44) 



(45) 



Bergman .epmputed the exact effective nonlinear conductivity of a binary medium in the dilute limit, for spherical 
inclusionga'Ea. His result takes the form 



Xexact = Xl + i (^2 - Xl){T - if 



^ d{d + 8) ^ ^ d{d~i) ^ ^ ci(3d^-10d+16) ^2 



d + 2 



d + 2 



3(d + 2) 



p + 0(p2). 



(46) 
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Compared to ( pq ) , expression (^) becomes exact when Xi = (nonlinear inclusion in a linear host) . It also becomes 
exact for d = 1, and in the limit d — > oo where the field is E = Eq in each phase so that = (x)- Moreover, it is 
exact up to order T^, which is consistent with its correct limiting weak-contrast behavior. For d = 2, it is exact up 
to order T^. However, the term of order is wrong for 2 < d < op. The reason for this misbehaviour is that the 
dilute limit requires an exact computation of (£'^}q (to linear order)Q, whereas we approximate it through a gaussian 
average. 

On the other hand, expression (E5) yields 



X^ 



Xi 



+ {iX2-Xi)iT- 1)^ + XiT^ [2(d + 2)-4T + T^]}p + 0(/ 



(47) 



once again a result less precise than ([45[). 

At the present time, since it requires an exact solution for the one-body problem, the test of the dilute limit appears 
to be the most challenging one for nonlinear effective-medium theories. 

This test is illustrated in Fig. 4 which displays x°> as calculated from (|3^) and (|2|), against the concentration p of 
medium 2, for moderate dielectric contrast and d = 2. The thick line segments at p = 0, 1 represent exact tangents 
obtained from (^6|). The tangent at p — 1 follows from substituting (l—p) for p and interchanging the indices 1 and 2 
in (^). Though not exact in the dilute limit, formula (33) yields tangents quite close to the exact ones. The marked 



inaccuracy of formula (|42|) near p — and p = 1 results in a lower value for x in the whole concentration range. 



D. Percolative behavior 

Before examining the predictions of Equ. ( |3^ ) near the percolation transition, I first briefly review the critical 
behaviour of WNL composites, and discuss a flaw of the Bruggeman formula in this context. 

An iefiulator/perfect conductor binary mixture undergoes a percolation transition for a critical metal fraction 
p = pc"- For WNL phases, the critical behaviour of the field fluctuations is now well understoodE3. In the limiting 
situation where £i <C £2, and for p < pc, one observes a behaviour 

e°cxel(p-p,)-^ x"fxxi(p-Pc)-('^+'''^ (48) 
on the contrary, for p > Pc, one has 

e"oc£2(Pc-p)*, X°«X2(Pc-p)('*-"\ (49) 

where s and t are the superconductivity and conductivity exponents, and where k and jp^MC the noise exponents 
which characterize an anomalous nonlinear susceptibility enhancement near the thresholcEJ'LJ. As a consequence of 
the exact inequality (£"*) > (-E^)^, the exponents k and k' are necessarily positive. 

In the Bruggeman formula, pc = 1/d, and £0 — £i/(l — dp) for p < 1/d, £0 ~ £2(1 — dp) /(I — d) for p > 1/d (£1 <C £2), 
whence s = t = 1. Reporting these expressions into (|3^) yields k = k' = 0, a result shared by all effective- medium 
formulae based on (p^). 

These values can be compared to exact bounds obtained for k and k' on the basis of the links-nodes-blob (LNB) 
motdel for electric percolation, now commonly accepted as a model for e.g. random resistor networks (RRNs). These 
ara2a 

{3d-4:)iy-2t + l<K<2{d-l)v-t, (50a) 
{4- d)iy -2s + l < k' <2iy - s, (50b) 

(for d < 6; the values for d > 6 are that of d = 6. Moreover t is not defined for d = 1), where ly is the correlation 
length exponent: ^ oc |p — Pc\~'^ ■ Though the latter information is absent from the Bruggeman theory, an inequality 
ly > must hold in any situation relevant to percolating systems. In (|5^), the lower bounds must not be greater than 
the upper ones, which leads to lower bounds for the exponents s and t themselves, namely 

s>(2-d)iy + l, (51) 
t> {d-2)iy+l. (52) 

In the Bruggeman theory, s = t = 1 so that only = - a problematic value by itself - is allowed, save for d = 2. 
Reporting v = Q into ( p0[ ) leads to k = k' = — 1, which is unacceptable. For d — 2, one obtains k = k' = 2;^ — 1. 
Compared to the values k = k' — of the above effective-medium theory, this gives v = 1/2, a non-conflicting value, 
though lower than the exact result 4/3. 
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Of course, the fact that Bruggeman's formula is a poor approximation near the percolation threshold is well-known. 
However, this discussion enlightens a fundamental inconsistency in the Bruggeman exponent values which goes beyond 
mere numerical inadequacy, as long as one wishes to estimate the properties of systems obeying the LNB scheme. 
It is to be noted that for d — 2, where no inconsistency appears, the exact equality s — t holds (a consequence of 
self-duality )e3 and is obeyed by the Bruggeman exponents. 

The above argument is to be brought together with another one by Bergman. Discussing the failure of a "non am- 
biguous" non-linear Bruggeman-type model, he concluded that bounds for fluctuations always refer to some particular 
type of microstructure, and that no such bounds exist which apply to all materialsu. 

Because of these problems, comparisons between effective-medium theories built on the Bruggeman formula, and 
simulations on RRNs are expected to be more significant in dimension d = 2 (another reason for preferring the 
two-dimensional case as a test-bench is that the bond percolation threshold pc = 1/2 on a square lattice is exactly 
reproduced by Bruggeman's formula). In Fig. 5 are displayed formulas ( ^3| ) and (^) against p, in a high contrast 
situation (e2/£i = 10000). The same trend as in Fig. 4 is observed: formula ( p3| ) yields-|a higher estimate. The values 
for £1^2 and xi,2 are the same as that used in Fig. 3 of a paper by Levy and BergmanE3, where simulations results on 
RRNs are reported, and compared to (|4^). The authors remarked that the height of the peak in at pc was badly 
underestimated by (^2[), which gives Xpcak — ^■'^ only. In their numerical simulations, for a 30 x 30 network, the peak 
height is found to be Xpoak — 0.86. In the infinite-size limit, this value is expected to increase somewhat. In Fig. 5, 
the height of the peak given by (|3^ ) is Xpoak — l-^' which thus compares fairly well to simulations. 



E. Other power-law nonlinearities 

In this section, potentials of the type (||) wherei 

<'(£;) = 2^^7+1 (53) 
7+1 

are considered. For simplicity, 7 is assumed to be constant in the material. The gaussian averages over the field in 
each phase are carried out numerically, with the help of a bi-variate integration routine. One integration variable is 
the field modulus E, and the other is the cosine E • E°/ {EE'^). In fig. ^ are shown plots for 

X° = {x{ElEy+^), (54) 

computed from gaussian averages, for a moderate dielectric contrast and various powers 7 > 1 ranging from 7 = 1.5 
to 7 = 5, against the volume fraction p of medium 2 (solid lines). An exponential enhancement of x" is observed with 
increasing 7. Its origin lies in the existence of a non-zero probability for E > E'^ in the composite, since E° = (E) 
by definition. These fluctuations are amplified by the nonlinearity, while that for which E < E^ are reduced. The 
larger the width of the probability distribution, the larger this enhancement, so that the peak culminates in the region 
p ~ Pc = 1/2. Once again, the estimates of the present theory are much higher than that predicted by the decoupling 
assumption 

xl = Y.P^X.{{ElEy)'^^^'y\ (55) 



(dashed lines), especially for large values of 7. This is understandable, since the functions to be averaged can be written 
w'f{E) = x{x)f{E'^), where f{z) = z^^^'^^^^/{l + 7). For 7 > 1, these functions are convex, so that {f{z)) > f{(z)). 
Therefore, the decoupling assumption always underestimates the fluctuations when applied to such potentials. An 
overestimation would instead take place with concave potentials. 



IV. CONCLUSION 



In this article, a theory for the nonlinear susceptibility of weakly nonlinear composites is proposed. This theory is 
based on a multimodal gaussian approximation for the overall probability distribution of the electric fleld, to linear 
order. The parameters which define this distribution (means and second moments, in each phase of the disordered 
medium), are obtained from Bruggeman's theory, whose limitations are discussed in this context. The present model 
provides for the first time an analytical estimate for the probability distribution of the electric field, for arbitrarily 
high dielectric contrast, in percolating media. The resulting effective nonlinear susceptibility is exact to second order 
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in the contrast, and close to the exact result in the dilute limit. Significant quantitative improvement is obtained 
on previous nonlinear effective-medium theories, even in the percolating regime, at least in the two-dimensional case. 
This study emphasizes the importance of accounting for the vector character of the electric field when averaging the 
local potentials. Improvements of the present theory could consist in finding a better approximation than the gaussian 
one for the components of the probability distribution of the field. More realistic distributions could be obtained by, 
e.g., extending the perturbative theory of Barthelemy and Orlandcj to high dielectric contrasts via self-consistent 
effective-medium approximations. 
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APPENDIX A: SELF-CONSISTENT ANISOTROPIC LINEAR THEORY WITH PERMANENT 

POLARIZATION 



With local potentials of the form (Uj), the homogeneized potential reads 



(Al) 



where 



9^J = 


fdn^ 


k-i kj 


J Sd 


















poi pO 





(A2) 

(A3) 
(A4) 
(A5) 



and the Bruggeman condition for e^^ is {fiij ) ~ Sij . 



Formula (Al) has been derived by means of a functional formalismEa, used with a trial potential of the formE2l 
W°{E) = (l/2)£^i£° i5j +P^Ei + c°. Rather than to give a demonstration which would comphcate the present article, 
the validity of (Al) is established by looking at its consequences. First, one can easily check that, for a binary medium, 
(Al) exactly reduces to formula (2.16) in the work by Ponte Castafieda and KailasancJ. Next, the average field per 
phase, {Ei)a = M°', is (using the symmetry of the tensors g and fig; /i itself is not necessarily symmetric) 



1 dWunjE^) _ . 

~ ^ — dP" — ~ ^'^^^ djk^Pk )■ 



(A6) 



From the definition of P?, one deduces that (M) = E^ , in agreement with the boundary conditions. Finally, the 
macroscopic consitutive relation derived from the effective potential, namely 



dE^ 



P^ 



(A7) 



is consistent with the equation 



{e^,,E,+P,) = {s%Mf + P^), 



(A8) 



with M" computed with (A6). This directly follows from the equivalent form of the Bruggeman equation, (sijfijk) 



^%{t^jk) = and from the identity /xg = % + fi'^k9ki{£u 
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FIG. 1. Field probability density function V\i{E\^/E°) 
= E ■ E°) in the Bruggeman approximation, for various 
volume fractions p2 of component 2. The dielectric constant 
ratio is 62 /£i ~ 5. Dots: p2 = 0.05 (violet, highest curve); 
P2 = 0.25 (blue,lowest curve). Solid line: p2 = 0.5 (green). 
Dashes: P2 = 0.75 (orange, lowest curve); P2 = 0.95 (red, 
highest curve). 
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FIG. 3. Field probability density function 7'||(-B||/S°) 
(i5|l = E ■ E°) in the Bruggeman approximation, for various 
volume fractions p2 of component 2. The dielectric constant 
ratio is 62 /ei = 1000. Solid: p2 = 0.12; dashes: p2 = 0.24; 
dots: p2 = 0.36. 
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FIG. 2. Field probability density function V^{E^/E°) 
{E± — ||E — iSiiE^II) in the Bruggeman approximation, for 
various volume fractions p2 of component 2. The dielectric 
constant ratio is 82/81 = 5. Dots: p2 = 0.05 (violet, highest 
curve); p2 — 0.25 (blue, lowest curve). Solid line: p2 — 0.5 
(green). Dashes: p2 — 0.75 (orange, lowest curve); p2 — 0.95 
(red, highest curve). 
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FIG. 4. Comparison between Eqs. (|3|) (solid line) and (^) 
(dashes). Effective nonlinear susceptibility vs. volume frac- 
tion p of component 2. The space dimension is d = 2. The 
parameters are: £1 = 1, £2 = 10, xi = 0.01, X2 = 0.1. The 
segments at p = and p = 1 represent exact slopes computed 
with expression 
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FIG. 5. Comparison between Eqs. (solid line) and ( p^ ) 
(dashes). Effective nonlinear susceptibility vs. volume frac- 
tion p of component 2. The space dimension is d = 2. The 
parameters are: ei = 1, £2 = 10*, xi = 10"'', X2 = 0.1. The 
segment at pj= 1 represents the exact slope computed with 
expression (|46|). 
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FIG. 6. Comparison between Eqs. (M) (solid line) and (|5^ 
(dashes) for power-law nonlinearities w" = xi^)^'^^^ / (7+1) 
Effective nonlinear susceptibility vs. volume fraction p of com- 
ponent 2. The space dimension is d = 2. The parameters are 
£1 = 1, £2 = 10, Xi ~ 0.01, X2 = 0.1. From bottom to top 
7 — 1.5 (blue), 2 (cyan), 3 (turquoise), 4 (green), 5 (black). 
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